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foreword 


The purpose of thi's manual is to ftoeument the' pjrpcedures and 
theory of oper^tjop flf the tAfjfDSAT trac|dng system at New England 
Diyision, U.3. Arniy Corp|i o'f En^inper§» Waltham, Mas§)achi|isetts. 

Thp manual is arranged generally by degree of detail, with 
ti\e sitjiplest'Qpe rating procedures fipst; instructions for normal 
^?i^rto-day operat^ori are- given in iSection I, whil^ Ip^onnation 
neecjed for prcjgram modification, file maintenance, and trouhle-^ 
shooting is in Sections II t VH and the jftppendiee?. All figures 
■peferp^d to in the text ^re ip Appendix B. 
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I. AUTOMATIC TRkONG SYSTEM OPERATION - OVERVIEW 

The Automatic Tracking System for receiving L’ANDSAT^'data at - 
New England ‘Division, Waltham, ' consists of a 15-foot dish antenna, 
a tracking pedestal, some pedestal control equipment,' and a' Data’ 
General NOVA minicomputer with various accessories. The--relationship 
of all these parts is shown 1n the subsystems diagram, Figure 1. 

Most day-to-day -operati on of the system will require very little 
action by operators, but full control of it and the handling of 
unusual situations require some knowledge of the programs and various 
information files that are kept on disc. Most operator action is 
taken at the computer terminal (Figure 7). (No card decks are needed, 
and the system is almost entirely separate from the IBM 1130 and 
Motorola equipment which are In the same room.) The operator may 
have to power up the NOVA computer (see section II) and start execu- 
tion of the programs which track the LANDSAT and store incoming data. 
Once the NOVA has the correct time of day and is executing the track- 
ing programs, it should do so continually until the operator interrupts 
it. These programs are cyclical, and If one of them is interrupted, 
it may be restarted later; that is, the operator may re-enter the 
cycle at one of several points. 
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The simplest procedure for tracking is as follows (ref^r to Figure 3) 

1. Power up the NOVA (see Section II). 

2. Set the fipVA's real time clock precisely (see Section III). 

3. Turn on all required control equipment (see Section V). 

4. Executp the program TRACl^ by typipg "TRACK" followed by a 
carriage return. 
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' •IT: ROUTTNE-'TO POWER UP THE'NOVA ' - 
SWITCHES- ARE, -UBELLED WITH R^D -JAPES. .(SEE FIGURE 7). 

1. TURN CONSOLE POWER SWITCH tO "ON ‘ 

2. TURN ON POWER SWITCH ON TEKTRONIX TERMINAL 

T ' t _ > 

3. TURN TELETYPE SWITCH TO' "LINE". 

4. TURN DISC POWER SWITCH TO' "ON" (clockwise). 

5. PUSH WHITE KEY ("POWER-ON")' ON MOViWG HEAD DISC' CABTNET' 

6. SET CONSOLE SWITCHES 0, 11, 12, 14, and 15 UP; ALL OTHER 
NUMBERED SWITCHES STAY ’boWNi 

7. WHEN GREEN LIGHT ON DISC CABINET COMES ON, LIFT "RESET" SWITCH 
AND THEN "PROGRAM LOAD" SWITCH ON THE CONSOLE. THE FOLLOWING 
DIALOGUE ENSUES: 


FILGH^E 7 mrSYS^ 

BDOS REV 3- 02 

bate ? 3 15 76 i. 

TIME CHiHiS> 7 13 15 0 1 

R - - r . 

CLEAR/ A/ V I' 

KR VSER) 

CLEAR/ A/V ) 
sirs, DR 
R 
R 

SfS* DR 
R '* 

NOTES: 

1. "J" means "RETURN". You type in the underlined characters. 

2. WHEN RDOS SYSTEM CRASHES, PUSH WHITE KEY ON DISC CABINET 
(OFF) AND GO TO STEP 5 ABOVE. 
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ROUTINE TO POWER DOWN THE NOVA (SEE FIGURE 7) 

1. TYPE RELEASE DP0 } ♦ ("0" IS A ZERO: "J” IS A RETURN). 

2. fUSH WHITE KEY ON DISC CABINET (WHITE LIGHT GOES OUT). 

3. TURN pFF FHD SWITCH. 

4. TURN DISC POWER SWITCH TQ OFF (counter clockwise), 

5. TURN TELETYPE' SWITCH FROM %INE^' TO "0FF% 

6. TURN TEKTRONIX TERMINAL SWITCH TO "OFF". 

NOTE: NOVA IS NORMALLY LEFT RUNNING ALL THE TIME. 

‘ 1 
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lU. SETTING NOVA'-s REAL TIME CLOCK 

f 

To track LANDSAT accuratley the NOVA's Real Time clock must - ' • 

be set to Coordinated Universal’ Time (UTC). Accuracy of one-fourth 
second -is sufficient. Two methods may be used; a manual one and - 
(eventually)' an automatic’ one.' 

' ‘ ’ 1 * ."J 

Automatic Method. Execute the program CL. Within 2-3 minutes the 
computer will signify completion by typing "R". If it doesn't Jt« 

means that it probably won't. In this case, use the manual method. 

- ■ '■} 1 

Manual Method. The NOVA provides for its clock to be .set by. the - 
teletype command STOP hh mm ss 

Where- hh^, mm, and ss stand for hour, minute, and second. 

Dial up the FTS number 8-323-4245 to get the National Bureau of 
Standards' audio time signal. When you have found out what time 
mark will be .coming soon (e.g., the next minute), use the STOP 
coimiand to prepare to enter that upcoming time, and hit CARRIAGE 
RETURN exactly when the time marker occurs. 

Before hanging up the telephone, you may check the NOVA's time .by 
executing the program PU which will send an audible pulse to the 
terminal every 15 seconds. 
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This annotated sampl^ of operator/computer d-j^log illustrates the 
method: 


* 


STOD 13 10 30 i 

r— r 


PUi 

IT'S NOW 13:10:34 
45 
60 
I NT 
R 

PUJ 

IT'S NOW 13;11:12 

15 

INT 


R 

LSI 1 

NEXT PASS AT 14:28:1{^ 


At this time, NBS signal is coming in by 
telephone. Return is pressed at exactly 13:10:30. 


Continue listening to NBS and- cpmpare teletype 
pulses tq telephone signal. 
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IV. HOW TO ENTER ORBITAL ELEMENTS INTO TRACKING SYSTEM 

xTo track' LANDSAT, -the system ;musti be able-to predict -when' the 

satellite will rise over the horizon and what azimuth and elevation 

angles to send to the tracking pedestal. To .predict those times 

and angles 'the' system is given a description of LANDSAT's orbit by 

means of the teletype or CRT terminal. This orbital information is 

contained in'the ^element set provided by the North American Air 

Defense Command, Ent AFB, Colorado.* The element set- comes via TWX 

twice a week and looks like the example in. Figure 5.- Eight of the 

elements in Figure 5 are important to our system. Their meanings, 
...» * * * . ’ 
formats, and units are as follows: 

1. EPOCH - - An arbitrarily chosen recent instant expressed 

i_ 1 

as a Julian date, at which the rest of this element set was deter- 
mi ned . 

XXX.XXXXXXXX (DAYS) 

2. NDOT0** - - First derivative of mean- notion + or - .XXXXXXXX 
(REVS/ DAY/ DAY) 


*Questions about NORDADC elements can be addressed to: 

SPACE DEFENSE CENTER 
(Cheyenne Mountain Office) 

ENT AFB, Colorado 

As of 9 December 1975, our contact person there was Capt, 
Tohlen, FTS 8-327-0111 635-8911, ask for ext. 3549 

** "0" stands for zero; "0" is the 15th character of the alphabet. 
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3. 10 - - Inclinatiori* XX.XXXX (DEGREE?) 

4. NODE0 - - Right Ascension of the Ascending Node, XX.XXXX 

(DEGREES) • ■ 

5. E0 - - Eccentricity. XXXXXXX (NO UNITS) 

Notice that the decimal point is not printed on NORAD message, 
but must be supplied tp system when you type it in. 

6. OMEG0 - - Argument of Perigee. XXX.XXXX (DEGREES) 

7. M0 - - Mean Anomaly. XXX.XXXX (DEGREES) 

8. N0 - - Mean Motion. XX. XXXXXXX X(REVS/DAY) • 

* 

The orbital element set is entered into the system by executing a 
program called "ELW", which stands for "^ement Writer", ELW 
is an interactive program which guides the operator in entering 
the numbers correctly. Because the numbers have many digits^ it 
is easy to mistype them on the teletype keyboard. Therefore, ELW 
echoes each number as It is entered and allows revision of that pne 
number- If no correction is needed, the operator types after "OK?" 
and enters thp next number, If a correction is needed, the opera- 
tor types and retypes the same number. An example of the opera- 
tor/computer dialog for thp element set of figure 5 is given in 
Figure 6. 
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V. POWERING' UP TRACKING EQUIPMENT 


- Power switches for the Data' General equipment,, the Sctentific/ 
Atlanta equipment and associated devices are shown on the photo-- ■ 
graphs in Figure 7. The order of turning switches ON is' as follows; 

Data General : 

' 1. Console power 

2. Disk power oh console 

3. Disk power on disk cabinet 

4. CRT' "I 

> reversible 

5. Teletypej 

6.. Decoder 

Sci entifi c/Atlanta : 

7. Receiver 

8. Synchro Display 

9‘. Servo Control ■- 

10. Digital Comparator 

’ in addition, the-' main power switch oh "the'^antenha pedestal 
concrete foundation must be ON,' and ahy''nnterl0ck''sw1tcheS' in'^the-'; 
pedestal itself must be closed.' . . ■ • 

fhere is one poWer switch on a plug stHp inside 'the S/A- cabinet 
of which you should be aware. 
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VI. IF SOMETHING FAILS ... 


Occasionally, a device will nialfimction and cause tracking to 
cease. Here are a few things to notice as you recover from the 
malfunction: 

1. Is the computer still up and running? 

If the lights are glowing softly, it is probably still run- 
ning and should respond to comnands from the terminal. If the 
console lights have stopped with some on and some off, or if no lights 
are on it has crashed. Go to Section II, Step 2, in the footnotes. 

2. Is the dish in the stow position (pointed straight up)? 

If not, tracking is either in progress or has ended abnormally. 
If it has ended abnormally, the command equipment must be returned 
to STANDBY mode. Go to next item. 

3. Are the two tiny (5/16") red lights (LED'S) on the Servo 
Control Unit lit? 

If so, the equipment is in PROGRAM mode. Put it in STANDBY, 
by typing “OFF" at the terminal followed by carriage return. The 
tiny red lights should go off. If that doesn't work, turn off the 
power switch on the DECODER (see Section IV) for 2 or 3 seconds. 

To put the antenna in the stow position, use the manual conmand 
unit, the cabinet immediately above the Servo Control Unit. Push 
the two square buttons on the manual command unit. The antenna 
should go to the position indicated on the two round dials on this 
unit. If it doesn't respond to the manual commands, something has 
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blown out - probably a fuse. Finally, .push the STANDBY buttons 
on the Servo Control Unit. 

4. If it looks like a fuse, has blown, zero in on the dif- 
ficulty by the equipment's behavior and appearance (e.g., lights 
out, movement in azimuth but not elevation, or vice versa) , lost 
power in Servo. Amps, etc.). Most commonly the center (30 amp) 
fuse on the concrete fbundation supply box is the one that is 
blown. 

Fuses have also blown in the Digital Synchro Display Unit' 
(back panel) and in the Servo Amplifier in the pedestal itself. 

5. The root cause of these fuse troubles seems to be back in 
the NOVA computer. When it sends bad data, the pedestal equipment 
gets overloaded. The NOVA runs into difficulties when RDOS is not 
functioning right, and recently it has appeared that RDOS gives 
problems when the files CLI.OL and TLOG are not "cleared". That 
is, their user counts in the system directory are greater than 
zero. This is remedied by the commands "CLEAR'TLOG" and "CLEAR 
CLI.OL" at the terminal. 

6. If the equipment all seems to work, but no signal comes in 
during a satellite pass, the system time may be set wrong. Also, 
check the system date. Perhaps the satellite has not been turned 
on by NASA at Goddard Space Flight Center. Usually, they turn it 
on by the time it reaches 10° elevation. 
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7. If the equipment has failed so badly that it can't be fixed 

by NED, the following service groups are available: 

Data General Corporation 
Field Service 
237 Riverview 
Waltham, MA .02154 
891-7024 

Tektronix, Inc. 

Field Service. 

482' Bedford" 

Lexington, MA 
861 -6800 

Sci enti fi c/Atl ant a 

Bud Lydon, Fred Leavett, or Dan Pioli 

Burlington, MA 

272-1256 
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VII. SYSTEM DESIGN AND OPERATION. IN-DEPTH VIEW 

The LANDSAT tracking system integrates a set of about twenty 
programs or subroutines (software), about ten disk data files, 
and several pieces of equipment (hardware). The inter-relation- 
.ships of the programs and data files can be seen in the flowchart 
in Figure 4. The hardware configuration is shown in Figure 1. 

In the flowchart, an information flow can be seen as well as a 
cycle of program executions. Essentially, the system predicts the 
satellite's position, tracks the satellite, stores and prints the 
data, returns; to the predicting program, and so forth. This 
cycle can be entered by the method given in Section I. However, 
from time to time, other operator action will be necessary. 

For example, the computer's real-time clock must be accurately 
set (see Section III). As of thi s' Wi ting, only a manual methods 
is available, and it has to be performed at least once a day for 
various reasons such as clock inaccuracy and system crashes. A 
better, automatic method of imputtirig time from a standard clock''- 
is being developed by the writer. 

The operator must' also inform the tracking-'system 'of the latest 
description of the LANDSAT' s orbit. This must be stored in a file 
called "ELEMENTS". The orbital information is contained in an 
eight-number set which is supplied to WCB under a standing- arrange- 
ment with the North American Air Defence Command (NORADC) in 
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Colorado. Twice a week, NORADC sends the element set to the TWX 
machine in Building 115S. A detailed description of how to enter 
the elements into the system is given in Section IV. 

•Normally, entry to the tracking system is by execution of the 
program TRACK, a FORTRAN program which calculates pair's of azimuth 
and elevation angles to LANDSAT from NED. TRACK starts with the 
current time and keeps, incrementing it until it calculates that 
the satellite would be above, the horizon.. In other words, it pro- 
jects into the future from the current time in, the computer's 
real-time clock. 

After the current date and time are input to TRACK, they are • 
converted to Sideral time. This is done by a subroutine called 
TCALC. . Sidereal timp is a relationship between the constellatipn 
ARIES (Y) and the Greenwich prime meridian. Specifically, it is 
an angle between the Greenwich .prime meridian and the inertial 
X-axis, which points toward the first point of ARIES (Escobal, p. 
20)*. This angle is denoted by 6. This angle is called the 
local sidereal time. Knowing the east longitude (Xe) of an ob- ^ 
server's .station and, Og (Greenwich Sidereal time) 0 can be easily 
determined. This is given by 0 = Og-- Ae, where 0< 0 — .2'TT 
(Escobal, p. 20, Eq. 1.26)'. 


*See Appendix F for literature cited. 
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To find the sidereal time, the Julian Date (J.D.) must be calculated^ 
The- Julian Date is a continuing count of each day elapsed si hce*- 
some- arbitrarily selected epoch. The epoch selected for LANDSAT 
orbital predictions in TCALC is January 1, 4713 B.C. Each Julian 
Date is measured from noon to noon; hence, it is an' integer 12 
hours after every midnight (Escobal, page'17). After the Julian 
Date and sidereal time are calculated, the next thing found is 
TSINCE, the number of days since the most recent NORAD EPOCH. 

TSINCE is then used to determine the unit vector pointing toward 
the satellite (see SGP of this text). These unit vectors are in 

i 

turn converted to azimuth, and elevation angles at the observer's 
station (degrees clockwise from north and degrees above the hori- 
zon, respectively) (see SRV of this text). These two values and 

the times at which they occur are then written on the disc under ' 

1 

* • M 

the file name PTAE. The Time of interest is then i ncrem ented' seconds 
or minutes by the routine INCT, and the. next values of AZ and EL 
are determined, etc. In this incremental fashion, the computer is 
able to predict the path of the satellite across the sky at the 
observer's site. ' , • 

The name PTAE stands for Paper Tape Azimuth' Elevation. The 
file can be transferred as it is to the paper tape punch by the 
teletype command, "XFER/A PTAE $TTP". This will cause a paper 
tape to be generated that is suitable for input to the ‘paper tape 
reader on the pedestal control equipment. 
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PTAE is a' disk file which comprises a time, a time increment, 
and many pairs of azimuth and elevation angles. The file.is 
ended 'With a special file terminator. An example of PTAE is shown 
in Appendix C. . 

TRACK calculates azimuths and elevations at TO-second incre- 
ments, so the angle pairs in PTAE are pointing angles for instants 
10 seconds apart. If these angles were fed to the tracking ped- 
estal, the antenna would jump quickly to the next position every 
10 seconds, the progress of the satellite is smoother than this 
jumpy motion, and it has been found that one-second incrementing 
is sufficiently small for constant satellite acquisition. There- 
fore, the program INTERPl is executed right after TRACK to inter- 

t 

polate ten angle pairs for every one pair In PTAE. Furthermore, 
INTERPl recodes the angles from ASCII characters to a binary 
coded decimal (BCD) format suitable for the electronic interface 
enroute to the conmtand equipment. The new angles are stored in a 
binary file called BCDAZEL, and the number of angles in'BCDAZEL is 
stored in the file MANGLES. 

The recoding is done by bit-mapping in the program INTERPl 
(q.v.); the assignment of angular values to bit positions is shown 
in Figure 8. A set of special interfaces, built by Robert Snyder 
of NASA Wallops, is used to route the BCD angles from the NOVA to 
the 1848 Digital Comparator (see Figure 1). Input to the 1848 is 
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in the form of Two 18-bit BCD words representing azimuth and 
elevation. However, the NOVA can output only 16-bit words, by 
way of the 4065 Digi'tal Interface. For this reason, it was , 
necessary to concatenate two pairs of 16-bit words into two'18-btt:? 
words as shown below (see also Figure 8).: 



18 -BIT SPECIAL 16 - BIT 

WORDS interfaces WORDS IN 
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■ .:TCAliC. • 


TCALC is a time handling routine which calculates the following 
three variables:' 

1. XJd - No. Julian days; it is used to find TSINCE And 
THETA (Escobal » Pgs. 20, 21, 22). 

2. TSINCE - No. of days since .most recent NORAD EPOCH. 

This is used by,SGP..ta find the unit vector pointing toward the 
satellite (ROOT). 

3. THETA - Sidereal time {measured in radians) - the 
angle between a line from the center of the earth to the first 
point of the constellation ARIES' (7C)s and the plane of observer's 
meridian,*' 



tHETft is used by SRV in determining Azimuth (A) and elevation (H). • 



EXPLANATION OF TCALC VARIABLES 


EP = 2442413.5 = Number of Julian days from an original EPOCH 
to January 1, 1975. This EPOCH is January 1, 4713 B.C. (days). 

TWOPI = 21T = 6.2831853072 (no units). 

LAMBDA E = East longitude from Greenwich to NED = 288.784332'* 
(degrees), ■( A© ) 

DTHDT = . 25068447 - Constant used *to account for one extra 
sidereal day for every tropical year (degrees/mim) . 

EPYR = DFLOAT (75) = An arbitrary year used as a reference 
(years). 

XJD = Number of Julian days (days). 

N = T(2)-l = Number of months in year up to last month (months). 

DAYS IN MO (1) = Number of days in each month (days). 

N = T(l)-1= Number of years up to last year (years). 

TSINCE = Number of Julian days at INSTANT, the time of interest 
(days). 

INSTANT = Future prediction times, or the times of interest 
(year, month, day,, hour, minute, second). 

DT = The number of. hours, minutes, seconds which T, the time 
of interest, is incremented for successive executions of TRACK. 

TU = Time since January 1, 4713 B.C.; used- to find THETA G 0 
(centuries). 
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THETA 60 = Greenwich Sidereal time at 0 hour of a particular 
date (degrees). 

THETA 6 = Greenwich Sidereal time (degrees). 

THETA -.Sidereal time at NED, Walthanij (degrees) . 

SAMPLE CALCULATION OF SIDEREAL- TIME 
August 23, 1975, at 10 hours, 15 minutes, 0 seconds; Number 
of hours, minutes, seconds expressed in minutes: DT = 6ls minutes. 

XJD = 2442648.5 DTHDT = .25068447 

TU = (XJD-241 5020)736525 = .7564271047227926 
THETA G0 = 99.6909833 + (36000.7689) (TU) + (.00038708) (TU)2 
= 33T. 648591601 5490 

THETA G = THETAG0 + (DT) (DTHDT) = 125.8195406515490“ 

THETA = (THETA G + LAMBDAE) = .9530184529430946 

radians 
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SGP 


SGP is a FORTRAN subroutine embodying a truncated simplified 

general perturbation theory for use in the determination of 

LANDSAT pointing elements. SGP computes osculating position, 

velocity and mean classical elements. SGP Is a first order 

analytical integration of the equations of motion including 

perturbations caused by the first two zonal harmonics of the 

geopotential. The zonal harmonic constants account for the 

effects of the non-circularity of the meridian cross sections 

of the earth. The. perturbations caused by these harmonics 

are independent of the longitude of the satellite. SGP is 

based on the orbital elements a, AXN, Ayn, 1» > and L 

which are well defined for all elliptic orbits except those 

that are nearly equatorial. For equatorial satellites, the 

elements' Axn and Ayn are ill -defined because of the indeterminacy 

* 

of the node, angle to which they are referred.- The SGP 
mathematical model is adequate to handle a majority of routine 
cataloguing. Accuracy is said to be better than one part in 
109 . 
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SRV 


SRV (Slant Range Vector) is a FORTRAN subroutine of TRACK 
which transfoYms the. orthogonal vectors and the time angle, 
THETA, from subroutine SGP into an azimuth/elevation coordin- 
ate system with the observer's station as the origin. Files 
of azimuth and elevation angles in this coordinate system 
describe the path of LANDSAT over a particular station during 
some interval. 
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TRACKING THE SATELLltE: PROGRAM LSI 

After TRACK has predicted the satellite's path across the sky 
and prepared a file of pointing angles, it chains automatically to 
the program LSI which will perform any of over six main functions. 
It is a complex multi-tasking- program which defies flowcharting, 
because program internal control shifts according to time as 
counted down by the Real Time Disc Operating System (RDOS) and 
according to real events in the outside world. 

Typically, LSI carries out the following main tasks: 

1. Schedules itself by looking at the smarting time of the 
upcoming pass. This time is the first number stored in the disk 
file BCDAZEL., 

2. Orients the antenna 1-1/2 minutes before the satellite 
rises. 

3. Starts repositioning the antenna second by second begin- 
ning at the instant the satellite rises; and simultaneously logs 
any data that arrives by way of the antenna/receiver/decoder path- 
way (see Fi.gure 1); and also simultaneously will accept corrections 
from the terminal to advance or retard some number of seconds. 

These corrections are made to improve antenna position. 

4. Restores the antenna to the stow (upright) position when 
the last angle pair in file BCDAZEL has been sent. 
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5. Dumps the field data that, have come in from core buffer 
to .a temporary disk file called "SDF" (Satellite Data File). 

6. Finally chains to a program called QD3 which will decode 

! ■ * 

I 

field data from binary to an octal format similar to one used by 
NASA at Goddard. 

Note that once TRACK and INTERPl have been executed for an 
upcoming pass, LSI can be run at any time up to one minute, 40 
seconds before satellite rise time. Execution of LSI after that 
causes problems which are signalled by a "W" being printed at the 
terminal,. One then has to quickly reset the system clock; execute 
LSI; and when the computer eventually types enter positive 
corrections that stand for numbers of seconds to enable LSI to 
catch up with real time. 

If further knowledge of LSI is desired, the program itself is 
the best source. The original source code is copiously annotated 
with explanations of individual steps. 
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OUTPUTTING DATA: QD3 and P3 


QD3 and P3 are programs that condition the raw field data received 

r 

by the ground receive antenna for disk storage or legible output, ' 

QD3 

Output to QD3 is the disk file "SDF“ which was produced by 
LSI immediately after the last satellite pass. Output from QD3 
goes to a temporary file "LS2DAT'* and a permanent file "STORAGE". 

These file formats are shown in Appendix C. 

The arrival time of each DCP message is recorded by LSI 
by storing a seconds counter with each message. This number of 
seconds is accumulated from the beginning of each pass, and QD3 
calculates a message arrival time by adding the number of elapsed 
seconds to the starting time. The arrival time (Y, M, D, H, M, S) 
is then stored with each message. The’ time used is Coordinated 
Universal Time. 

P3 

Legible output of DCS data is obtained by executing P3. Input’ 
to P3 is from temporary disk file "LS2DAT". Essentially, the 
program examines each message for the platform ID number. . looks 
up the ID in a table, and decides how to interpret the data on the 
basis of indices in the table. These indices then direct program 
control to appropriate subroutines for calculating decimal numbers 
and attaching labels. The kinds of parameters handled by P3 are 
shown in the sitelist. Figure 9. 
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APPENDIX A 


HISTORY AND BACKGROUND OF LANDSAT PROGRAM AT NED 

Since the Industrial Revolution In the 1800 ‘s, the rivers of New 
England have been developed to supply water for power and trans- 
portation. As new means of transportation became more economical, 
both railroad and highway systems were built along the banks of 
the rivers to service .the expanding needs of the industrial, com- 
mercial and urban centers. Structures, such as buildings, roads, 
bridges and dams have restricted floodways to such an extent 
that considerable property and environmental damages have occurred 
during moderate and major floods. Notable floods of November 1927, 
March 1936, September 1938 and August 1955 have demonstrated the 
need for flood control to prevent these natural catastrophes. 

At the direction of Congress, the U.S. Army Corps of Engineers 
developed a comprehensive plan of protection for each river basin 
after a careful analysis of all water resources. Protect! ve works 
generally consist of a combination of channel improvements, dij^es 
and/or floodwalls at major damage centers augmented by upstream 
flood control reservoirs. Many of these reservoirs contain addi- 
tional storage reserved for other uses such as water supply, con- 
servation and recreation. The Corps has built 35 flood control 
reservoirs, 37 local protection projects and four hurricane barriers 
in New England at a total investment of some $300 million. 
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To achieve optimum operating benefits from this comprehensive 
protection system, the New England Division requires hydrologic ‘ 
data such as river, reservoir and tidal levels, wind velocity and 
direction, barometric pressure and precipitation. 

In the past this data was collected from field observation and 
relayed via telephone or voice radio. It took several hours to 
compile and assess the data in this manner. With the need for 
timely and reliable information increasing, the Corps began develop- 
ment of new methods of data collection. 

In 1970, the Automatic Hydrologic Radio Reporting Network was 
placed in operation. This ground-based radio relay system con- 
sists of 41 remote reporting stations, and a central control at 
Division Headquarters In Waltham, Massachusetts. This network, • 
under computer programmed control, collects and analyzes, in real 
time mode, information which is essential for flood regulation.. 

The remote reporting stations are strategically located in five 
major river basins and at key coastal points, with each contribut- 

ing to a detailed, comprehensive hydrologic picture. 

LANDSAT 

In June 1972, NASA entered into a contract with the Corps for an 
experiment to study the feasibility of using the Earth Resources 
Technology Satellite (ERTS or LANDSAT) for collection environmental 





data fv"om Data Collection Platforms. (DCP's) which are installed 
at 27 locations throughout New England. Many are situated at 
existing U.S. Geological Survey gaging stations. 

Since July 1972, LANDSAT has been relaying river stage, precipi- 
tation, and water quality data from DCP's via the Goddard Space 
Flight Center to the U.S, Army Corps of Engineers, New England 
Division, in near real time. This is the first resources satel- 
lite designed to' obtain data from the planet Earth exclusively for 
planning, design, operations and research of land and water re- 
sources. 

THE NED GROUND RECEIVE STATION 

Since any operational satellite configuration serving an urgent 
function like flood control should include ground receiving stations 
at all major user locales, NED, with NASA support, constructed and 
-is now operating an inexpensive semiautomatic and easily maintained 
ground receive station as a follow-up to its original study. The 
Division is now able to receive hydro-meteorological data from data 
collection platforms in the field directly at its headquarters in 
Waltham, Massachusetts with no time delays. The software to drive 
the antenna system has been developed with the intention that the 
antenna operate in an unattended mode .automatically over nights and 
during weekends and holidays, with a computer controlling all 
processes. 
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The major objective of the program has been to compare the effec- 
tiveness of the LANDSAT-Data Collection System (DCS) with existing 
systems in aiding our watershed management functions. 

Data collection platforms tested by the Corps have performed suc- 
cessfully in all seasons including the winter months and also dur- 
ing significant flood events, transmitting near real time operation- 
ally useful data- for our flood fighting missions. 

The satellite proved invaluable in April and early May of 1973 
and 1974, monitoring flooding in Maine Rivers. LANDSAT relayed 
data from five river points in that State to aid the- New England 
Division in the coordination of the flood emergencies. 

The successful testing of the LANDSAT Data Collection System at 
the New. England Division should encourage serious consideration 
of the institution of an operational satellite data relay system 
on a Corps-wide basis. Such a system appears to be more ^cost- 
effective than conventional ground-based data relay. 

The .New England Division is also making a study of satellite 
imagery to determine its usefulness in planning, designing and 
managing water resource systems. To obtain an overall broad cover- 
age of ground conditions, imagery studies and measurements are- being 
made of fluctuations in river, lake, and reservoir stages as well 
as tidal changes, icing of water surfaces, location and depth of 
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snow cover, moisture content of the soil, and'water quality para- 
meters. 

FLOOD CONTROL OPERATIONS 

Data received at the New England Division's Reservoir Control 
Center from either the Automatic Hydrologic Radio Reporting Net- 
work or the LANDSAT Data Collection System is compiled by computer. 
This is augmented by information from other sources such as the 
National Weather Service Meteorologic and River Forecast Offices 
and the U.S.' Geological Survey. Experienced engineers and hydrolo- 
gists at the Reservoir Control Center analyze the data for timely 
operation of dams and hurricane barriers, and then issue instruc- 
tions to operating field personnel. 

Flood Control reservoirs, local protection projects and hurricane 
barriers built by the Corps in New England have been responsible 
for prevention of almost $300 mill ion in flood and storm damage. 
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FIGURE 2. LANDSAT ANTENNA AND PEDESTAL 
DCP ANTENNA IS SHOWN TO THE RIGHT. 
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WNOS<»T*a - DCP INFORMATION SHEET 
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15 APRIL :9T6 
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STATION NAME 


7147 ST. JOHN RIUER AT NINEMILE BRIDOE, ME. 
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7201 PEMIGEUASSET RIVER AT PLYMOUTH, N.H. 
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APPENDIX D - GLOSSARY 


ADR - Aniog to digital recorder. Typically a Fisher-Porter or Leupold- 
Stevens recorder, equipped with a telekit. 

Azimuth - Horizontal angle measured clockwise from north. 

BCD - Binary .Coded Decimal. 

Chain - in programming, a call from one program to execute another, 
thereby terminating its own execution. 

Coordinated Universal Time - an observer's local mean solar time plus 
the number of time zones the observer is west of Greenwich observa- 
tory, corrected for aberrations in the spin of the Earth. 

Crash - (v.i.) - to cease functioning. Syn. bomb. 

CRT - Cathode Ray Tube - Specifically, the Tektronix 4014 terminal 
connected to the NOVA. 

DCP - Data Collection Platform - Field installation used for sensing - 
parameters, encoding data, and transmitting data to satellite. 

Disc (or disk) - medium for storage of data in the Data General Com- 
puter. Refers to twenty-surface disc pack and drive which is a 
peripheral device to the computer. Files on the disk are divided 
into two directories, "DP0" and “USER". Most system programs are 
in DP0, and most user programs are in USER. 

j 

Elevation - angle above, the plane of the observer's horizon. 

Flowchart - diagram that shows flow of control in a computer program. 
Elements shown are input, output, initializations, processes, 
decisions, and connectors. 

Julian Date - an arbitrary benchmark that is a continuing count of 
each day elapsed since some particular epoch. 

Multi-tasking - in a computer several program tasks competing for 
devices and the central processor on a priority or queued basis. 

Octal - refers to a number system that has 8 as a base. 
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Sidereal Time - the relationship between an observer's meridian and 
some inertial coordinate system, for example, one based on the 
constellation ARIES, 

Real Time Clock - device in the Data General Nova Computer that con- 
sists of a crystal controlled clock and associated DG system soft- 
ware that are used (1) to keep track of date and time of day and 
C2) to provide for low resolution timing. 

Tracking - keeping the antenna pointed at the satellite, and in 
conjunction with that, logging any incoming data. 

Universal Time - see Coordinated Universal Time. 
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RLW S/K 16/C LSI T1 la T3 T5 WG AHT GPRB/L TASICOALL FMTaB 
fORT.U;DCtETE/C GARB 
R ' 

TYPE LSI 

C LSI CALLED BV INTEPPl OR EXEC'D BV ITSELF 
C UAITS FOR. SATELLITE RISE TlflE - ZMIN, 

C ORIENTS ANTENNA^ TRACKS SATELLITE^ LOGS BATA 
C ACCEPTS CORRECTIONS^ BUMPS DATA AT END OF, 

C PASS, STO»J5 ANTENNA IN UPRIGHT POSITION 
C TURNS OH AND OFF CIRCUITRY IN COtWAND 
C EOUIPMENT, CHAINS TO 0D3 
C TDB 14 NOU 75 

COMPILER DOUBLE PRECISION 

-DIMENSION IDf3^IT(3),IT2(ll),IT3(ll ),IT5fll^ 
COMMON/KBLK/ITKll) 

EXTERNAL TUTa,T3,TS 
COHNON/KEY/KEYl , KEVa , KE V3 , KEV4 

C K£Y-"S APE USED TO PASS MESSAGES AS FOLLOWS: 

C KEY! T1 TO LSI 

C KEV2 T1 TO Ua 

C UG BACK TO T1 

C KEV3 TS TO T1 

C KEY4 T1 TO UG TINE COUNTER 


6 


10 


OOMHOH '^lELK/DAVSiriMOMa ) 

DATA 0AVSiriH0/3l.,a8.,31*,30.,3l.^3O.,3l*,3U,3^>. ,31., 
♦30. ,31./ 

ZM2N-1. S/1440. 

CALL DFILU("SDF'‘,IER) 

CALL CFILUi-SDF*,a,lER) 

CALL OPEMtS^-BCDAZELM^IERj 
CALL APPENDa5,"TL0G-,3,IER» 

CALI 0PEMUa,*5TARTANGLEM,IER > 

CALL 0PEm7,»MAMGLESM,IER) 

CURDATE^O. 

CORTiriE*0- 

KEVa-0 

KEV4»0 

CALL TIHE<1T,IER) 

IF<IER.KE.1)TVPE “TINERRMER 

CURTinE*DFLOAT<n( 1 ) H^DFLOATdTCB ) 1/60. DFLOAT( IT( 3 ) ) 
♦ 3600. > 

CALL DAT£(ID,IEP) 

IF(IEP.NE*UTYPE *DEMER 

ii-ioa>~i 


DO'S 

CUPDTE • CURDTE^DA VS INMO ( J ) 

CURDTE-CURDTE^DFLOATdD^,a> ) 

1F< < IDO )/4<4.EG.ID( 3> I.AND. ( IDU ).GT.a) )CURDTE«CURDTE^-1 
CllR0TE»CURDtE*CURTlNE/24. i CURRENT JUL TIME SINCE I JAN 
READ! 5, lO^END-100 )TS1NCE,DT 
F0RnAT<F13.9,F6.1> 


. 1F(CURDTE.GT,364.)TYPE PRESET EP AND EPVR IN TCAtC OH 
t4AN, 1 ♦SEE WOOD PAGE FOR MORE INFO" 

IF(CUR0TE.LE.(TSINCE-2NIN))G0 TO 50 >TA$K SCHEDULING 
TYPE *U- ;T00 LATE FOR CURRENT PASS 

CALL CHAIN! 'TRACK •SUMER) 

IF<J£R^HE.X)TVPE •LStRDSlER 
CALL EXIT 


C 

C 


AHCLES^Ta ORIENTS, T3 GATfCRS DATA 
TS ACCEPTS CORRECTIONS 


50 READ(7,51)ITl<a) ;M0, TINES TO EXECUTE 

•: NOTE THAT ITUa) IS THE LOC THAT IS MODIFIED BV TS 

51 FORMAT 116 J 

XTIME- (TSinCE-IDINT aSIHCE j j r24 , 

ITK4)-XTIME tST^^TING HOUR 

ITl<5>-C<TlME-DFL0AT{m(4>M*36O0. ;SECOHD UITHIN HOUR 

ITll6)-3 

ITt(7^«DT 

ITKin-lOO 

CALL FOTASK(DUM,Tl,ITl,lER,-n 
IFdER.NE.l >TVPE TQl^IEP 
IMIH»m(5)^G0 
I5EC*m(5)-II1IN*60+.l 
URnEa0,60)ITl<4),IMIH,ISEC 
URITEt IS, GO I m f 4 >, IMIM, ISEC 
60 FORtlATC" NEXT PASS AT - 12,* ;M2^ : M2 ) 

lTar2*»l ;OPIENT ANTENNA ONCE 
XTIME- i i TSINCE-ZMIH ILaNT(TSIIKE-2t1XM ^ UB4 . 
IT2C41«IDINT(XTIME> 

IT2^5 !•» XTIME-DFLOATf IT214 ) > )f 3600. 

na^6i-3 

na»7 wo 

naai)«200 

CALL F0TASKiDUN.TadT2,!EP,-l t 

TF< lER. HE. UTVPE “FOaMER 

IT3!2)-i 

lT3MJdTit4i 

IT315 W1T1<5 J 

IT316 W0 

IT3(7)*0 

IT3(11)-300 

CALL FaTASK(DUM,T3,lT3.2£R,-l j 

IFUER.HE.l JTVPE 'F03 MER 

lTS(a)-l 

IT5M)*ITI 4i 

ITSl5)‘ITl«5l 

IT5(S)-10 

IT5(7>-0 

ITSai)*SOO 

CALL FQTASK<DUI1,T5,IT5,IER,-n 
IF (lER. HE. 1 iTVPE ’'F05MEP 
10NE«0 

CALL REC^KEYurONE'^ 

CALI FDELVUetW ;UAIT FOP ANT TO PCM STOU POS. 

II-l 

CALL ANTdl, 11,11,11, II> ;TURN OFF CIRCUITRY 
OALL CLOSERS, lER/ 

CALL CLOSED 7, IER> 

CALL CL05Eda,IER) 

WRITE aO, 70) 

WRITE! 15,70 « 

70 F»>HnHT<,' » 

CALL CHAINrOD3.SVMER) 

IF< lER. HE. I >TYPE •CHER*,IER 
100 CALL EXIT 

CALL CLOSE (IS, lER) 

END 



TVP£ T1 

C TDB 9 DEC 75 
TftSK Tl 

COfinON/'KEV/KEVl , tCEVa> KEV3 , KEV4 
copinoN/<Buc/iTiai) 

DIfENSION 

ie*9 

KEY4*KCV4+1 
IF<KEV3.GE.0)GO TO 10 
KEV3-KEV341 
ITl<e)»ITt(21*l 
GO TO 25 

10 READ BIH«?VCS)l<t> 

IFtra).EO.-l)00 TO so ?EMD of FILE7 
READ BINARVlS)a<J),J-a,'4) jNO, GET 3 MOPE 
CALI AWT(ia).H2),I(3),K4M0) jSEMD THEM TO ANT 
1F^KEV3.LE.0)GO TO 85 jADJUST IF NECESSARV 
tCEy3*»CEY3-l 
m(a)-iTi<8)-i 
GO TO le 
as CALL KILL 

30 CALL XI1T<KEVa.l,$50J jTELL UG TO DUMP TO M5C 

lOHE-0 

CALL REC^FEV2,I0NE> jUAIT FOP UG TO FINISH 

CALL XMTiKEVl,t,8S0) jTELL LSI EOF HAS BEEN REACHED 

CALL KILL 

SO TVPE “XPITERRI* 

CALL EXIT 
END 
R 

TYPE T2 
C TDB 8^7/75 

TASK T2 
OINEMSION J(4) 

CALL APPEND(15,*TL0G*,?.I£R) 
ie«o 

^ READ BlNARY(ie)J ; GET FIRST ANGLE PAIR 

C FROM ’STAPTANGLE* 

IPOINT*! 


GO TO 100 ;G0 SEND PAIR 

^"<5^ ANTEHMA to GET THERE 
DO 80 1-1,21 ;SEHD ANTENNA CCU TWICE 
READ BlNARViiau 
IPOINT-2 
GO TO 100 
CALL FDELYCai) 

CONTINUE 
GO TO ai 


100 

lao 

21 



URiTEas,iao> 

URITEaO,120> 

FORNATi^ ,',2> 

GO TO <10,l5.25,3S),IPOlHT 

iPOirrT-3 
GO TO 100 
CALL FBELVCai) 

CONTINUE 

READ BINARVCSJJ iFWST LOOK ANGLE. 


IPOlNT-4 
GO TO 100 


jFWST LOOK angle, 

THEM CONTROL GOES TO Tt 


35 URITE(10,I30) 

WRITE(15,130' 

130 FORPf^TC- 

CALL KILL 
CALL CtOSEaS.ICR) 

END 

R 

r/PE T3 

TASK T3 
CALL 

CALL KILL 
END 
R 

- TVPE HNT 
•TlTL ANT 
.EHT ANT 

4EXTD •CPVL^.FRET 
*HREL 

1—167 ;L0 AZ 

J-I-H ;HI AZ 

K-I+a ;LC> EL 

L«I+3 ♦HlEb 

lM+4 

FS.-5 •PROG/ETDEY SUITCH 
FS* 

ANTiJSR e^CPVL 
STA 3 RTN 
IDA 0 «M 3 

NOU 0 0 S2R ;TUPN OFF PEE€STAL? 
JMP OFF ;VE5 

LDA 0 ATE *H0^ SEND ANGLES 

STA 0 TTB 

IDA a CIH 

NOR 5 LDA e 3 

LDA 1 «TTB 

COh 0 0 

AND 1 0 

con 0 0 

I5Z TTB 

DOAS 0 WC 
INC 3 3 

INC a a szR 

JMP MOR 
IKt LDh 3 PTN 
JSR e.FRET 

off: LDA 0 BO 
DOAS 0 DUC 
JMP Bk‘ 

iDATA AREA™ - 

DUC-4S 

ATBJ.+l 

B51 13S777 lOBl-^OBS 
133777 ,OB140S4 

B3i 1E7777 ^OBl40B3 
B5i U7777 jOB140B£ 
bo: lB04iB2 

RTN:,-. 

CN4t -4 
TTBi.-. 

> r 
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TYPE UG 
.TITL UG 
•ENT UG 

• EXTN .UIEX, •REC, •IXMT, .TASK, .AKILL* 
.EXTD •CPVL,.FR£T 

.conn KEY 4 ^ LABELLED COMMON AS IN L5 

.Txm 1 

.NREL 

sDEUICE CONTROL TABLE 

IDDCTtDCT 
DCT # ♦ " • 

1B7 

ISR43 

j INTERRUPT SERVICE ROUTINE 
;ReSPOMl« TO ONLY ONE INTERRUPT 
;FROn 4065 INTFC, tOADS S UORDS IN ft 
;ROU AS FAST AS DECODER PROVIDES THEM. 
;TIfMMG MATCHES DECODER'S E5^ACTLY 

ISR43t NIOC Due 
STfl a URTN2 
STft 3 URTN3 
LDA 0 SVNC 
JSR e.TOBUF 

LDA 1 CNS 




nORJ DIA 0 Dve 
JSP 9.T0BUF 
JSR TtIP 
INC 1 I S2P 
JMP nOR 
JMP OUT 

.TOBUFJ TOBUF 
TOBUF: UDA 2 «.MBP 
STA 0 0 2 
INC 2 2 
STA 2 e.MEP 
JMP 0 3 


TMRS LDA 2 CM 15 
INC 2 2 S2R 
JMP ,-l 
JMP 0 3 


^'gl'DA^0^«BnESS jfiET TIME COUNTER • 

f!$ i j STORE TIME UITH MSO 

INC 3 3 
$TA 2 «,nBP 


SUB 1 i 
LDA 2 URTN2 
LDA 3 URTN3 


NIOS DUO 
•UIEX 


umcf.-. 




SYHCt 12214 
CNSs -5. 

CNIS: -IS. 

CN3lA0i-2l00. 

CTRi 

.ntPJ MBP 

UGU5R t.CPYL 
STA a Aca 
STA 3 AC3 

SUB e 0 :GEN A 0 
LDA 2 CN3i 60 
STA 2 CTR 
LDA 2 PBUf 

STA 0 0 £ ;INIT BUFFER TO ALL O'JS 
IHC 2 £ 

ISZ CT» 

JtIP .-S 

LDA 0 DUCN ;DEFINE 4065 DIGITmL I 0 BOARD TO SYSTEM 

LDA 1 IDKT 

.SVSTM 

.IDEF 

JMP (fc.ERT 

.SV5TM 

-GDAV 

JPP «.EPT 

2 3 ;NE/T 7 UHU STuPE ■/P.MO.DA'-' IH SDF 
LDA a PBUF 
STA 2 riEP 
STh I 0 2 
INC 2 2 
STh 0 0 2 
IHC 2 2 
STA 2 MBP 
.5V5TH 
.OTOD 
JtIP 0.EPT 

HO‘l 2 3 ^riEXT 7 LINES STOFE HR, WIN, SEC IN SDF 
LDA 2 0.P1EP 
STA 3 0 2 
INC a a 
STA 1 0 £ 

IHC 2 2 
STA 0 0 3 
IHC a 2 
STh 2 f.MBP 
LDA 0 AMESS 
SUB 1 1 
HIOS DUC 

.REC jUAIT HERE FOP LhST LOOK AtIGLE TO BE SENT 

LDa 0 DUCN 

.SVSTH 

. iRnu 

JftP 0.EPT 

LDA 0 ASDF i POINT TO FILE HAftE 

SUB I 1 ;DEUICE CHARS 

.SVSTH 

.APPEND 6 

JtIP 0.ERT 

LDA'0 PBUF 

nOVZL 0 0 jBYTE POINTER 
LDa 1 C4S00 ;BVTE COUNT 
.SVSTM 
»URS 6 ; 


DUMP THE BUFFER 





J«P i.ERT 
.SVSTf! 

.CLOSE 6 
jnP «*ERT 
LIM 2 AC2 
LDA 3 AC3 
LDA e AnESS 
SUBZL 1 i 
.xnr 

JKP «.CRT 
JSR «.FRET 

iDATA AREA 

AC2J 

AC3j,-, 

AIIESSt.GADD KEV,1 ;POINTS TO 2ND ELEMENT IN LABELLED COrtHON,*' 
BMESSt .GADD K£V,3 ;POIHT$ TO TIME COUNTER >SECS^ IN T1 
ASDrt .+1*2 

•TXT ^USERtSDF* ; SATELLITE DATA OUTPUT FILE 
C4200r 4200. 


DVCNME 

DVCM2 

riBPJO ;P10UADIE BUFFER POINTER 

.ERT:ERT 

ERT: .SVSTM 

.ERTN 

PBUF:BUF 

BUT? *BLK 2100. ; BUFFER FOR 3O0 NSGUTUDS/flSG 

.END UG 


R 

TYPE TASKCALL 
.TITL TASKCALL 


.ENT .TasK,.:<MT,.>^MTM, *REC^.|i:iLL..OUE>!,.TOULD 
♦ENT .OUREL, .OUKIL* .QTSK^ .PRI, .SUSP, .PEND» .IDST 
.ENT .TIDS, .TIDR. .TIDK, .TIDP, .AKILL, .HSUSP 
.ENT .ARDV, .APEND^.AUMPD 

• > XMTT, XNTTL) , RECC > KILL , TQUEX, TOULD . TO^RL^ TOi >KL 
.EXTN T0T5K,TPRI,TPEND^TIDST.STID>RTID.KTID,TIDP 

♦ EXTN TAOL^TAPEN.TWJNP 



TPRI 


•SUSP - 

J$R 

e. 

TREND 
.IDST - 

JSP 


TID5T 
♦TIDS ■ 

JSR 

«. 

STID 
-TIDR • 

JSP 

e. 

RTID 

.tidl: - 

JSR 


KTID 
♦TIDP • 

JSP 


.hKILL 

• JSR 


TAKIL 

.ASUSP 

- JSR 

A 

■TAPEN 
.ARDY - 

JSR 



TAUNP 

;EOUIUALEMT RDOS CALLS 

.PEND » .SUSP 
.APEND . .ASUSP 
♦AUNPD ^ .ARDY 

.END 

R 


GTOD 
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IDRIGINAI] PAGE IS 
OP POOR 9,UALE133, 



fOM TRACK OTl TCALC SGP SW IHCT SEfll EXftHn PEAK - 
GARB/L EORT.LB;DELETE/C GARB 
R 

TVP£ TRACK 
C TBB 1DEC7S 

COnPILER DOUBLE PRECISION 
DIHENSIOH TT(6).1D<3) 

CALL FGTin(I,J,KJ 
TT(4)-I 
TT<5)-J 
TT{G)-0 

CALL DATE(ID,IER) — 

TTa)«X0f3) 

TT(2>»I0{1 ) 

TT<3)*1D<2) 

CALL OTHTT) 

END 

R ‘ 

TVPE TCALC 

C TAR R/RR/7R 

CONPILEP DOUBLE PRECISION 

SUBROUTINE TCALC(T,TSINCE, THETA, LAMBDA E) 

REAL LAMBDAE 
DIMENSION T( 6 ) 

COMMON 'IBlK'DAYSlMMOdE/ 

EP-2442778.5 ;l JAM 76 
TUOPl -6. 8831853072 
DTHDT».350E3447 }DTH£TA/DT 
EPYR*DFL0hT(76> 

XJD*EP+(nn-EPVR>t36S. jADD 365 DA-'YF THRU LAST VR 

DO 50 ;DAYS IM P10NTH THRU LhST NO 

XJD*>^JD+DhySItlMO(U 
50 • CONTINUE 

N-Tt 1>-1* ;N-L*hST VR 

DO 100 I»?5,N jCHECk: FOR LEAP VP5 THRU Lh5T YR 

lF(I/4J:4TrtE*I iGO TO 106 

XJD-XJD+l. 

100 CONTINUE 


N-Tri> ,MOU LOOk. AT THIS VR * 

IF<N/4<4.ME.rOGO TO 200 :IS THIS MOT h LEAP 
1F(T(2).LE*2*^G0 TO 200 ;ARE ,UE &EVOMD 2/25- 
XJD-XJD+U ;VES, ADD A LEAP DAV 
200 XJD-XJD+T(3) 

C JULIAN DATE AT INSTANT 

TSIHCE-XJ0+T(4)/£4^-MT(S>+T(6 V 6 O 4 J/1440. -EF 
DT-T<4)1:60. + T<5 ) 4 -T <6 )/« 0 . ;(HR,HirKSECS » AS WINS 
TU- C XJD-24 15020 . 0 > / 36S2S * ‘ 

THETA G0-D«OD( (99.6909833-f< 36000 ,76S9tTU J+.0003370?f 
•TU«2i,360.) 

THETA G«DNOD( ( THETA GO+*r>TmMDT ) ,360 • > 

C SIDEREAL TINE IN RADIANS 

THETA- (DNODCf THETA G+LAMBDA E^360 ♦ ) )T7UOPI/360* 
RETtRH 
END 
R 

TYPE SENI 


C TDf 2/28/75 

CONPILER DOUBLE PRECISION 

double precision function SEnrcEE,XNN,XII) 

CONPUTES THE WEAN (KOZAI) SEHI-NAJOR AXIS OF A SATELLITE 
YY-. 3333333333 
X42- *00108248 


XNU- 11467. 25238 
AA-(Xmj/XHNI:T2)ttV7 

DD— 1 * 5rXJ2T U I . /AA J/( ( DSOPT( I . -EEJl^tZ n tt3 ) 

DD«DDTa.-1.5*(DSINail > 

SEfll ' AA^ ( 1 . + YY I^DD- YY U>D^ T2 ; 

RETURN 

END 

R 

TYPE ITKT 
0 TDB 2 28. 75 

COMPILER DOUBLE PPEOISIOM 
SUBROUTINE INCTIT,DT^ 

DinENSiort.Tce*' 

COmON/IBLX/DAVSINfiOaa) 

, IFfDT.GE.eO. JGO TO 600 

T(6)»T(6;+DT •UiCR SECONDS 
IF<Ti6).LT.60. 3G0 TO 350 
TC6i*T(6)-60. ;RESET SECONDS 
600 IFfDT.C^.60. )T'5)-T^5)+DT^60. 

IFfDT*GE.60. iGO TO TOO 
'T*^5>-r^5)+l ilNCR F11NUTES 
700 IF^T^5».LT.60. iGO TO 350 

T(5’i.nS)-'60. ;R£SET MINUTES 
T(4^rTi4)fl. ;INCR HOUR 
IFtT(4.MT.24, GO TO 350 
TM**T(4i-24 iPESET HPS 
:«712) ;PTP TO NO 
Tl3)-T^3'^-l. ^rNCP DAY 
IVR-Ti 1 t 
I LEAP '0 

IFa.EQ.a.hND.IVP 4*4. EG. ISP iILEAP-1 
DAVSINM0f2'-2e+rLEHP 
IFiTi'3j.LE.DAVSiriNCM I I iGO TO 350 
T‘3»-X. j RESET DAYS 
DHVSIMNOiE »«2S. ;PESET FEB 
T*'2i-T«2^+1.‘ jlncR MO 
IFa»2J.LE.12. «G0 TO 350 
T^BJ^l. ; RESET MO 
T 1 jiriCR VR 

350 return 
END 
F 

TYPE EXHNN 
r TDB 2^28 75 

COMPILER DOUBLE PPECISIOH 
DOUBLE PRECISION FUNCTION EXANnarW,ECC) 
CONFUTES ECCENTRIC ANOMALY USING KEPLER'S EQUATION 
TWOPl -6. 283 1853072 
EXANM-DMOr^XMn.ruOPl ) 

DO 10 I-1,S0 
«A*ECC*DS1N<EKANM> 

DELn-XMM-EXAN^I+AA 
ZZ- 1 . -ECCtDCCS ( EX ANN ) 

DELE’-DELM 7 224( r .SfDELfH/ZZUAH) 
IF(DABSIDEL£>-1. >30,30,20 
20 DELE-DELE/DABS<D£LE J 

30 EXAHN»EXAH#HDELE 

I F »' DABS ( DCLE )-*O0OO0n40,10,10 
10 CONTINUE 

40 CONTINUE 

RETURN 
END 
R 

GTOD 





TYPE on 

C TUB 1DEC7S 

COHPILER DOUBl£ PRECISION 
. SUBROUTINE OTlfT) 

OIRENSION r<6],TS(6> 

PERL 10, 1(1,N0,HDOT0,rj0,«OOe0,KODEf1,L«, Nrt,MDOT«. LLONG, 
tL0,NODDT,J5,J3.nU,IS,NOD£S,timSDA E 
EXTERNAL SEIU 
INTEGER VR 

COIWON EPOCH, YR,HO,HODE0,OnEG0,NDOT0, 
*pn,£«.ln.NOII£W,OHEGn,L«,W1,MDOTin,E0,H0,r0,L0,A0, 
tELOHG, LLONG.EXLMG, OBEGL.TRUEU, RHAG. ROOT, NODOT, ONGOT. 
#UX,OV,UZ,PX,RY,RZ,R!>OTX,RDOTV,RDOT2 
ccmmyiBLK/Mysimx 12 > 

«30*^3l 

c6Hy'.0174S3232Sl j DEGREES TO PADS 
1COUNT«0 
J2*. 00 108248 
Q— 1. 

JQ’l 

AE»1. 

IMAHV-0 

X2*0, 

2 « 0 . 

X*0 


IW*1 

TUOPI*6.283i353072 
CALL APPENDUS, 'TLO<5‘, j,IER) 
IFflER.NE,! ;TVPE ’TLtRMER 
CALL 0FEH(5,-ELErlENTS!,l,ieR) 
CALL IfILLK 'PTAEMER) 

CALL CPIUN'PTAE'.2,IER) 

Call 0P£N(4,*PTAE',3,IER) 

C It«T 

oneSScoE?™- 

lO-IOlCONO 
IIODEO'NCiDEOtCOHV 
ti0*H0»coNu 
N0"t«0»TUOPI ’ 
HDOrO-NDOTOUTgOPI 
CALL CLOSE 15 .IEP) 

IPdER-NE,! iTS'PE *CE*,IER 
LO«W04NODEO4OMEGO • 


A0*SEI1I<E0,nO,I0) 

iS5ft^^*l^!jS'*.Jf*^*'E-'^A0*<l.-E0**a# . >txz 

HODOT»-TEnP.tDCOS ( TO ) 

2 . -2 . St < DSIN( 1 0 ) )»*2 J 


SDT-ie. 


DT'BDT 


XI -DT 


e 

69 


CALL 0PEN(8,*C00RDM,IER) 

IPtlER.NE.t JTVPE *COE*,IER 

LWdI ,SNPHI, CSPHI , G2SHPHI 

LAUBOA E LAMBDA C-. 89077 ;E,L, CORR TO COINCIDE UITH 


iPASS»t norad prediction 

URITE( 10,69) 

URITE<1S,6S) 

DURATION' 


3S Tsa)»Ta) 

aa CALL TCAUCiT,T5INCE,THETA,LA«BDrt E) 

^ TSIMCE*TSIHCE-EPOCH ;HO,^OF DAYS SINCE HOST RECENT NORAD 

CALL SCP(TSINCE) 

IF(X2!Nr .«2CSPHI, SNPHI,CSPHI,02SWPHI ) 

lGO-0* 

Z*0 

X>0 

IMAt«Y* IHAN't’+I 
J0>1 

^ 

STEL-90.00 ;STOU ANGLE 
0RrTE<4,339<H,STEU 
IFdriANY.NE.IPHSS'&O TO 990 
CALL CLOSE* 4, lEP/ 

CALL CL0SEt8,IER) 

CALL CL0SE<15,IER) 

CALL CHAIN dltlTERPl. S' <‘,IEP. 

IF' lER.NE.l )T'/PE -LScE" 

CfiLL EXIT 

CHECK TO SEE If SATELLITE IS ABO'JE THE H0RI30M 

339 IF(H.LT.O. WO TO c5A 

IF(G>.G£.e..0P.IFLA'3.E*I.l )G0 70 190 
IFLAG‘l 
DO 150 1«1,6 
ISO T»ri*T5(I' 

DTxSDT 
H-Q 
X2*DT 
GO TO 20 
180 Q'H 

lOUER- I 

CHLL PEHK( IHAN'f', H,2,y, ICOLP'T , A, DT, tC-"£R > 

IFLAG»0 

IFtJO.EO.l <yRITEaF,i3F.-Tt4),T<5),T(S) 
IFiJO.ECi.nyPITE(10,l?SiTt4),Tf5),T(S» 

135 FORNATax,2(F3.0'* ‘ ),F3.0,2) 

J0»3 

TSS-TSlNCEttPOCn 
IF*IU.E0.n«RITEl4,333)TS5,DT 
333 F0P.HAT<IX,F14.?' 1>,F6.1 ) 

URirEl4,32?)A.H 

339 FORf3AT(2H A,FS,2,1HE,F6.2, •: ') 

IU<2 

290 DO 300 I«l,6 

300 TStfTd) 

CALL IHCKT.DT) 

CO TO 20 

250 IF(Q.GT.0. )DT-BDT 

lF(X2.eO.0)GO 70 887 
IFiO,GT.0)X2-DT 
387 0-M 

tU>l 

GO TO 290 
EHD 
R 

CTOD 
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TYPE SGP 

C BY TDB — REU a 14 75 AT 1430 
COWILER DOUBLE PRECISION 

c thiI'^r^tim^conputIs^tellite position using a simplified 

C GENERAL PERTURBATIONS HETHOD, CLASSICAL MEAN ELEMENTS APE 
C IWUT* AND POSITION, UELOCITV, B OSCULATING ELEMENTS APE 
C RETAINED 


REAL 10, Id ,N0, MDOTO, MO, HODEO, NODEM, LW, NN, NDOTM, LLONG, 
»NDOTe,L0,NODDT,ja,J3,nU,IS,NODES 
EXTERNAL EXANM 
INTEGER VR 


- INPUT PAPAMETEPS 

COMMON EPOCH, YR, no, NODeO,Ot1€uO,NDOTO, 

#A«,EM, I«,NODEn,Ofl£GM, LM,NM, NDOTM,EO,NO, 10,L0,A0, 
*ELONG,LLONC,E:<LNG,OI«';L,TRUEU,RnAG,RDOT,NODDT,OI1GDT, 
#UX , UY , UZ, PX , RY, R2 , RDOTX , RDOTV , RD0T2 
ja-. 00103248 
J3»~.O0OOO2SSa 
HDOT6-0. 

• AE-1. 

■ MU-1, 

TU0PI*6.aS31853073 


COMPUTE TIME UARIAHT MEAN ELEMENTS AT TSIHCE 
TT^TSINCE jTIME SINCE, EPOCH lOAVS- 

DM»HO17T4t)DOT0*TTt»a*NDOT6i:TTtl:3 jCHG III MEhM hNC'WAL. . 
DOMEG'OtlGDTITT ;D ARG PER 
DNODE»NODDT»TT jD hSC MODE 

LM-DMODi >L0+DNtOOMEC,«-DNODE »,TUOPI 1 ;I1EAN ORBITAL LHGT'JDE 
OMEGH'DriODi m.'ME&0+DOMEG1,TUOPI i ;A.P. 
NODEH-DMOD<tNoDE0+DHODEj,TUOPI 1 ;RA OF AN 
!M«IO ; INCLINATION UNCHNGD 
SINI’DSIH' IM' 

COSI«'DCOSaM> 

Nri»IW+2.tND0T0*TT+3.4!ND0T6*TT U2 

AM'AOi: ( ( I HO/HM UT. 3333333333 m2 ) 

Ef5*l.-AO'ANtU,-E0) 

IFl£M)ie, 10,20 
10 • EM -0.00001 


COMPUTE AND APPLY LONG PERIODIC TERMS ISUESRCPTD' “L* • 

30 TENPL ■ iJ3/ jaiKAE/AnjISINI/d.-Emia) 

AXHL-EM*DC05(0MEGM) 

AVNL-EMTDS IN I OMEGN . 5TTEMPL 
ELONG- DSORT.i AXNL**a+AVNLT« ) 

• • 0MEGL-DI10DC(DATAHa<.AVNL.AXNL ••sTUOPDj PRESERUE OUAD 
C LONG PERIODIC ON I- ISt 

' , LLONG- DMOD < < LH- , 25 tTEMPLTAKNLT 1 3 . +5 . TCOS I >/ < 1 . +COS Hi, 
OTUOPI), 

C SOLUE KEPLER'S EOUATIOH AND OTHER TUO-BODV FORMULAE 

C LONG PERIODIC ECC ANOMi 

EXLNG>EXANM ( LLONG-OMEGL-NODEM. ELONG / 

C TRUE ARC OF LATITUDE! 

^ TRUEU-2, »DATAM(DSaRT( < 1 .BELONG ) '< 1 .-ELONG n| { DSINf . 5* 
«EXLKQ>/OCOSC .StEXLNGO > ItOMEGL 

RflAG-AfWd.-ELOHGtDCOSCEXLNGM ;R SUB L 
C TRRdSUERSE COMPOHEMT OF UEl UECTOR 


f* o o •.“> r* c* 


RUDT-DSORTvMUJAflTCl .-ELONG**a ; >»' I./RMAG) 

C PADIAL COMP OF "EL UECTOP : 

PMODT-DSORTIMUFAH *tELOHG.'RMAGtDSINfEXLNG) 

COMPUTE AND APPLY SHORT PERIODIC TERMS 

TEMPS-. 254J2«« AE^‘, Hl1t( l.-ELCNGt tad HT2 {PERTURBATION 
SIHaU«DSIN<2.tTRUEU> .CONSTANT 

C052U-DCOSi2.tTRUEU) 

RHHG-RIIAG+TErPSiSimttatCOSeUnAPIti l.-EL0NG4*2> ) 
TRUEU*DMCD< dRUEU- .S*TEMPSt> 6. -7. tSIHITt 2 >tSINau ',TUOPI > 
1S-IM+3.TTEMPS»SINIIC0SI»C0S3U 
. NODES-HODEM+S.nEMPSICOSIlSINau 

COMPUTE. OUAHTITIES POP OUTPUT 
SNODE-DSIIM MODES ' 

CNODE-DCOS< NODES J 
SINI'DSIIUISj 
.COSI*DCOS>1Sj 
SINU-DS lN'TPUEUi 
C05U-DC0S«TRUEU* 

. UNIT 'lECTOR POINTING TOUhRD ShTELLITEJ 

SEE P.104 IN ESCOBAL, ’METHODS OF OPPIT DETEPMINATION* , 

TO CHECK UAUUE6 

UX-COSUtCNODE-SINU4SNOD£*COSl 

.Uy-COSUI'SNODE+SINU4CriODEJtCOSl 

uz’SiNutsim 

V‘ — SIllUtCNODE-.OSUISNODEtCOSI 
W»-SIHUI5NODE4COSUtCHODEICOSI • 

U2-C0SU4S1NI 

UX-SIHUSHODE 

UY— SINItCNODE 

:J2-:0SI 

P\-RMAGTUX 

RV-PMAGNJV 

PZ-RMA'JtUZ 

C PPOTX-RMGDT lUX+R' 'DT N.'M 

C • RDOTV-RIIGDTIU'.’+PODTtUV 
C RD0T2-RMGDTIUZ+RHDTHC 

C RD0T»D50PT'RD0TM*2+RDCir't*2+PD0T2W2 ' 

PETUPN 
• END 
P 



2 - 2 " ^ly 


TVPF SRU 

c SRV (SLANT RAMGEv UECTOff ) TDB (fi/aT/TS) 

C COfWENTS AFTER LINES IN THIS SUBROUTINE ■ 

C ARE EQUATIONS IN APPENDIX OF ESOpEAL, 

•NETHODS OF ORBIT DETERMINATION* 

COMPILER DOUBLE PRECISION ■ 

SUBROUTINE 'SRU(THETA,H,A,GlCSPHI,GacSPHI.SHPHI,CSPHI, 

^^lEALS^e, IM,Ne,NDOT0,M0,NODE0,NOPEM,Ltt,NN.NDOTri,LLONG, 
«L0,NODDT,LX,LV,tZ,LXH,LVH,LZH 

INTEGER VB 

COMMON EPOCH, VR,M0,NODE0,OMEG0,NDOT0, 

BAM, EM. IM, MODEM , OMEGM , LM , HM, NDOTM , EO , HO, 10, L0 , AO. 
tELOHG,.LLONG , EXLNG, OMEGL*, TRUEU , RMAG.RDOT , NODDT ,0«GDT , 

*UX. U V , UZ , PX , R V , RZ , RDOTX , PDOT Y ; RD0T2 
■ TUOP I *6 -283 18530 73 
X*(GlCSPHI>lDC05(TfCTAf' {1A.62 
‘ Y*(G1C5PHI)*DSIN(TICTA))1A.63 
RHOX*PX+X ;1A.6S 
RHOY-RV+Y ;1A.63 
RH02-P2+&aSNPHI ;1A.7A 

RHOH-DSOPT(RHOX**3+PHOYT*24RH02*B3j jIh.?! 

L.-'»RHOX-'RHOH ;UHIT mecTOR FROM SITE JA.T3 TA SATELLITE 
LY-RHOV/RHOH ;DITT0 1A.73 
-LZ«PHOZ RHOH; DITTO 1A.7-) 

COSTH-DCOSiTHETh; 

SINTH-r>SIH< THETHJ 

LXH*LXTSHPHnC0STH+LYtSNPHU5INTH-LZ*0SPHI ;IH.7S 
IVH • -L:'.3 S IMTHt LV T C OSTH ;D ITTO 

L:H»LXT(OSTH»CSPHI*LVtSIMTHI.OSPHr+LZTSNPHI ; DITTO 
cOSH-DSORT«’i .-LZHTT3' 

H*DATAH'LZH- COSH • ;IH.T6 IN MOCiI' 

CORPECtlOM FOP PEFPhCTIOH FOLLOUS ;C0URTESV OF RALPH pHfS (iSFC 
H*h+.00O?«DCuS'H; -'CDSIMi H )+DSORT' .04*i DSlrliH • 'M2' > 
H»3b0.*H.'TU0PI ;♦ OP - DEGREES FROM HORIZON 
A’DHTAlta- LVH.-LXH ' jlH.77 
A»3S0.TA TUOPI ; DEGREES CU FROM HOPTH 
IFJA.UT.O. 'A-A+3S0. ;nDJUST COOPDIHhTE SYSTEMS 
RETURN 
END 
R 
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TVPE PEAK 
C JHB 

COMPILER DOUBLE PRECISION 

9UBR0UT I N£ PEAK ( IMAf IV , N, Z , , 1 C OU^T , A , DT , lOUER J 

C^ILL APPEMDil5.“TLrjG'*,3,IEP- 

ICOUNT-lOOUHT+1 

IP(H,LT.ZJ'iO TO, 333 

Z-H 

RETURN 

333 IFCv.EOa^GO TO 1 
m-tn^^Nv+i 
NU(1-ICOUNT-l 
ZMUri-NUM 

ZNUflWZNUntStDT • 60* 

rniouER.EO.0 jgo to lo ^ 

irfZ*EO*0 .OR*X.EQ*0HJRITE< 

10 ir^z,EO*0 •OR.^.EO.o 'ijRiTEf i0,a£2aM,z,:<nun 

FOPHAT‘ tt>,^n,S:<^F7*2,2X,F5.0/ 'I 
X-i 

IF«Z7t:Y7rs"*HJ^ir.H,LT*0^C>^LL FCHrtfK‘’UHlT20*5U‘' i 
1 CONTiriUE 

C^^LL CLOSE 1 IS, lER' 

RETURN 

END 

P 

CTOD 

3. 23 '76 16520:9 

F 



TVP€ n 



\ 

m 

I 

KO 


21 

16 

4 

Id 

5 
18 

19 

11 


2d 

3d 

4d 

45 

48 

Gl 

Gi 

62 

C 

54 

63 

C 


RHC 8 JAN 76 

DINENSION IDX< 14,25) , ( 6 ) > ITEMPSf 12 > 

C0Wt0N/JBU/NUP1(24),SUM(12),X<S),JBIN<9),LBIN<3^ 
COWOH/ JBLK/ J HAHE (12J,CMK 

CO«fK)H/J8LK/COND,DOX,TEhP,FH 
CALL DFILUMDATA'.IER) 

* CALL CFILW<*DATAV2,IER) 

CALL APPEND(S,*DATA-,3,IER) 

CALL APPEND(6,“TL0CiS3,IER) 

CALL 0PEM<7,*INDXM,1ER) 

CALL 0PEN(8,-DIRC*,1,IER) 

CALL 0P£N(9,*lS2DAT%l,IERi 
IFaER.rC.l JTYPE-OE- MER 
WRITE (6, 21 J 
URlTEae,21) 

FORriATUX,-PID -,2X,'STATI0N NAHE^9X,*DAT£%6X, •EST\ 
IX-d 

READ(7,16^IDX 

F0RI1AT(I4aX,Il,U,12A2) 

IFLG-0 

READ(8,ld,END«15d)IPD 
FORMAT (14) 

READ ( 9 > 1 2 , END* 11)1 TEMP 1 , IP I D , ITEMP2 
F0RMAT^6A2,I4,12A2) 

IFaPID.EQtIPD)L)RlTE(5,19)2TEMPl,IPID,ITEnP2 

iraPID,£O.IPD)IFLG-IRlG+l 

F0Rr)ATaX,6A2>I4,i2A2) 

GO TO 6 
PEUIND 9 

IF(IFLO,£0*e)G0 TO 39 
CALL CLOSEfS.IER) 

C ALL OPEm S , - DATA M ,IER ) 

W 15 IJM,1FLG 

READC5,2d,END-145)IV,If10,IDD,2HH,NM,;SS,IEPR,IPlD,WM 

FORN AT fIl,5I2,Al,l4,24in 

IFdERP.EO. 17696)00 TO 15 

CALL UTII ( lY, mo, IDD, ISS ) 

00 30 1*1,29 

, rraDX(i,i),EoaPiD)Go to 40 

CONTINUE 


IND-IDXC2,1) 

DO 45 J-3,12 
K-J“2 

NAME<K)*IDXfJ,I^ 

WRlTEi 6,48 )IPID,NAME, INO, IDD, iV, 

URITE< 10,48 )IPrD,NAHE,IMO,IDD,lY,IHH,Iiri 
FWmTax,14,lX, 12A2,IX, 12, ^ ^M2, VM2, IX, la, • ! Mi 

70(61,62, 63, 64,65,66,67, 68,69 )IHD 
yRITE(6,51 ) 

VRITEaO,51> 

F0Rf1AT<3x,*TE5T GET') 

GO TO 15 
CALL ERTDA 

IF<CWC*E04l,)G0 TO 125 
UR1T£(6,S4 )5 
URITE (19,54)5 
FORWAT C 3X, * STG- • , F6. 2 > 

GO TO 15 
CALL ERTDA 

ZF<CHK«EQa)QO TO 125 


URITE(6,55 )S 
WRITE (10,95 >S 

55 FOPMAT ox, “PRC- ■ , F6 . 2 ) 

GO TO 15 
64 COND-0 

DOX-0 


TEMP-0 

PH-0 

CALL UTOLV 

IF(CHK.EQ.l. iGO TO 125 
URITE^ 6,96 )COHD,DC»X, TEMP, PH 
URITE a 0 , 96 )COHD , DOX, TEMP , PH 

56 FORRATOX, "CD- • ,F6a,2X, 'DO- " , F60, 
?43X,*UT--,F6*2,2X,'PH-SF6*3) 

GO TO 15 

69 CALL SHOP(HUM, END, DEPTH) 

) C IF(CHIC*EO.nGO TO 129 

URITEaO,5TJDEPTH 
URITE<6,57 iDEPTH 

57 FORMAT < 3 X , - UE S - • , F7 , 3 > 

GO TO 19 

66 URITE I'b, 58 J 
URITE< 10,98) 

58 FORMATOX^-CRREL" ) 

GO TO IS 

67 CALL ERTDA- 

. CALL SNOP(nUM,IND*DEPTH» 

C IF.(CHK.E0.1)G0 TO 125 

URITE(6,59)5,DEPTH 
URITE (10,99)5, DEPTH 

59 format ( 3X , " STG * ” , F6 . a , 2X , “ IJES , F 7 . 3 ' 

GO TO IS 

68 CALL ERTDA 
COND-0 
DOX-0 
TEMP-0 
PH-0 

CALL UTOLV 

C IF(CHK-E0,1 >G0 TO 125 

URITEC6,60 >S,OOND,DOX,TEnP,PH 
URITE(10,6O iS,C0ND,D0X,TEf1P,PH 

60 F0RMATQX,**STG-%F6,2,^,43X,'CD-*,F6.l,2X,‘D0-',F6,3, 
7/,43.<,'UT-‘’,F6.2,2X,‘‘PH-'*,Fe.3> 

GO TO IS 

69 riRlTE<^6,52) 

UPITEa0,52> 

52 F0R(1AT<3X, ’'DEMO. 5ET* > 

) GO TO 15 

125 Wn£‘.6,130^ 

MPITEaO,130^ 

130 FORMATt 6X, " IHMALID* ) 

19 CONTINUE 

39 DO 152 IQ- 1,29 

2F(lDXil,ICl>*E0.1PD>G0 TO 154 
15£ OOfiTlNUE 

GOTO 4 

154 DO 156 JQ-3,12 

KO-JQ-2 

196 NAMEUCO )-lDXt JO, IQ > 

IFaFLG.EO,0)URlTE(6,16O)IPD,NAME 
IFCIFLa. EO . 0 JURITE ( Id, 160 )IPD, HANE 
160 F0RMATUX,I4,1X,12AS,4X,*N0 REPOPT') 



Ol-3<^ 


HS REUlMD S 

GO TO A 

150 CONTIhUC 

CALL CLOSED S>IERJ 
CAU CF X LUC "DATA MER^ 

CALL CFlLUC'DATA-^a^lERj 
CALL CLOSECe^IER) 

CALL CL0SE(7,IER) 

CALL CLOSE (8, I ER; 

CALL CL0SE(9,IER) 

URIT^C 10^2001 

aee forhat </////) 

CALL CHAirK • TRACK. 5C^".rEP> 

IF(IER.NE.t)TVPE*C£- MER 
STOP 
END 
R 

type ertda 

SUBROUTINE EPTDA 

C0nn0NMBLK/f<UM(24 ),LBItK3 v JElNi 9^XlS ‘,bUMU2 1 

COfinON/ J ELK / J.K.S,F,JX.iy. n A(1E U £ > . UHK 

COnnOM/ JBLK/COMD . DOX • TEMP, PH 

JJO 

5*0 


F-1. 


30 


10 


32 


31 

25 


2& 

20 


40 


Jl-1 

J2*Jl+2 

K«-3 

T-0.005 
DO 10 

IFtNUMf I J.OiT.? ^ CO TO 40 
K«K*3 

CALI DCBPvi LEIN. hum n.jj) 

DO 10 L-1,3 

J-L4K 

JBIfKJ I'LPIHCL r 
M-2 

DO ao 1-2,9 
IFa-6i31, 32,31 

f^io.kf 

T-0.005 

Kl’6 

TFtJBirUI Jl20*EB,20 

Ti-F^a.tca-Ki ) 

T-T+Tl 

IP<F:U0.-T)4O,26,26 

S-S+Tl 

CONTWJE 

Jl-4 

F-.ei 

IF(J2M >30^40, 40 

RETURN 

CHK-1. 

RETURN 

END 


R 

TYPE BXHAl 
C D C 1/1/7S 

subroutine BtNAl 

C0J^/JBU/J,K,S,F,JX,IX,NAHEC1^ 

connote 

W) 10 I«l,3 


IX-a^JX-O-D'fG 
ir(NUfiaX>-7.M5,45,55 
45 K-K+3 

CAU DCLPViLElH.KUm r/),3) 

DO 10 L-1,3 
J-L+K 

JBIHtJ i-LBlM(L> 

10 continue 

55 RETURN 

END 
R 

TYPE BINEl 

SUBROUTINE BIMEl 

C0NN0N'JELK"NUinf24\SUin( 12 ^>U3 JEIN« 9 ^LEIth 3> 

COmON^JBLK/J.K,S,F,JX,rv,NAf1EU2),CHK 

C OtmON' JBLK • C Of ID , DOX ^ TEflP, PH 

DO 45 1*2,9 

n*i-t 

1F( JBttUl ) 140,10,40 
10 IF (1-6)20,30,30 

20 XfMi-aUfS-I) 

GO TO 45 

30 X^NV.5U^9-I) 

00 TO 45 
40 

45 CONTINUE 

ri*o 

Du ZO N-1.2 
N*2UX-ca-Ni 

sunrn »-o 

DO 60 L-l,4 
LX*41H“(4“L) 
suntn'*suM(i^^+*-{U< » 

X 1 LX > ■ 0 

60 CONTIHUE 

IF(SUNCN’-10/80,7O.70 
70 L-e 

30 CONTINUE 

PETUPN 
END 
P 

TYPE PCBPV 
C D C i '1 '75 

SUBROUTINE DOBRV- LBIILNUM,JJ i 
DINEHSION LBIN.4J 
DO 10 1*1, JJ 
10 LEINaj*0 

IF(NCTn5,5,35 
35 DO 30 KK-l,NUtt 

' DO 15 IM.JJ 

IFtLBINU M20,2O,15 
IS CONTINUE 

20 DO 25 

25 LBINUI-O 

LBIN(J ^■l 
30 CONTINUE 

5 return 

END 
R 



ORIGINAI^ PAGCIS E-1 1 

OF POOR QOALSai, 


TVPE UTOUV 

SUBHOUTIHE UTOLV ^ 

COnOON/aBLK/'NUn(34 ), sum IE ),X< 8 ) , JFIN<9 1 , L8IH( 3 ) 
COWlOH/ JBU/ J . K, S . F , JX, I X , NAME (12 ; , CMK 
COnnOH/4BLK/CONO,DOX,TEnP,PH • 

DO 25 JX*1,6 
CALL BlNAl 
CALL BINBl 

as coHTiNue 

rF(tc-6n«.30,3e 
3e COHD-0 

D0X«0 
PH*e 
TEMB-0 

COHD* 1 . 54: ( IM *sun ( 4 ) +1 04SUn ( K tSUM ( 2 ) ) 
DOX-0.02T(10«4Slfl1(U )+10»SUI1(ia^*SUri{31i 
PH*0.014*(l«0»SUn(S)-HC«SUH(6)+SUH(7l J 
TEfiP*0.04»ao0<sun(8Hi0*suri(9j+suna0)' 

RETURN 
10 CMIC-1. 

RETURN 

END 


p 

r/P£ UTII 

C RHC 13 JAN 76 

C COtlUERT:? ZULU TIME TO EASTERN STANDARD TINE (EST) 
SUBROUTINE UTII « I V, IMO, It'D, ISS > 

INTEGER DA'/SINHO 
COMHOri. IBLK-'DAYSr»rtO(ia) 

DATA DAVSINNO'31,29, 31, 30, 31. ?0, 31, 31,33,31,30,31-' 

iFass,GE. 30 )Mfi«nn+i 

IF'Mn.GE.eO'IHH-^IHHtl 


rF(MM.OE.60;Pl[1*Mn-60 

IFaHH.LT.0,IDD«IDD-l 

IF -s IHH . LT . 0 1 1 HH» IHH+a^ 

IFllDD.LT.l <tllO-IttO-l 

iFafio.LT.i)ino-ino»i2 

IF 1 3 DD . LT • 1 .• I DD -DAYSINHOf ItIO ,* 

IV«lV+70 

RETURN 

END 


R 

TYPE SWOP 

SUBROUTINE SNOPrNUPI, I,DEPTH-» 
DIflEMSION NUM<24) 

CRCT»6.9S 

IF(I.E0.S)CRCT*3.92 

D6»64tNUrK 16)4Stt«i«( 17 )+NUmi8I 

D7-64tNUHtl9>+8«NUM<a0)+NW1<ai ) 

DEPTH- (0?-'D6 >455 . 36-CRCT 

RETURN 

END 


TVPE.LSTGE 

C SUBROUTINE TO CALCULATE DISCHARGES AND UALIDITY FRON 
C UALUES OF STAGE Ar*D STATION KUNBEP 


SUBROUTINE LSTGE ( I SH , STAGE . 0 , NAME . DA, J ) 

DINEHSION DlSCH(4AJ.NAMEa0) 

COHNON.'KBLF./ LS 1 US ) 

DATA LS 1/7 1 47, 7101 ,7230.7071 , 7272,7356 ,7201, 7233, 
$6063, 7106,7304,7345,7254,7335,7506 ,6504/ 


DO 85 N-1,16 
1F(LS1(N).E0.1SH/G0 TO 8 
85 CONTINUE 

GO TO 285 
■S NN-N-1 

CALL FSEEKU3,HHi . 

C T-INITIAL stage in TABLE} 0-STAGE IHCPEHEHT IN 2HD LINE- 
REAPU3, 3 'MANE , DA,T, 0, DISCH 

3 FORMATreX, 10A2,24 a,F6.0,12X,F6.O,6:{'(3X,9F8.OU 

C J'O FOR OAUID OR J-1 FOR NON "hUD SThGE 
J*0 

IFt 0)50,50,21 ;F0P I-UP1KY FILE, IHCRNNT (01 WILL BE A 

21 IFtSTAGE-r;5l,2E,23 

51 0— .005 
GO TO 285 

22 DO 60 H-1,43 
S-STAGE -T-NTCl 
R-CS+Q'/O 
Z-DISCHCN) 

P-DISCHlfN-ll-Z 
IF'P >24,60,65 

65 IFlS 170,70,60 

C INTERPOLATE ON PhTINC TABLE 

70 IFi 2 152,52.700 

52 1F>R)51.152,152 

152 0-P4-P**l.5' 

GO TO 50 

700 Q*Z+P»P+,005 

GO TO 50 
60 COHTINUE 

•: extrapolate c-n fating table 

24 0*Z+P»(2-DISCH(N-n 14.005 

2S5 J-1 

50 RETURN 

END 
P 



2 L- 3 -S '8 


TVPE 0D3 


C 

C 

c 

c 

c 


TDB le rWRCH 76 ' . . 

PROCRAM TO COMUERT LANDSAT MSGS TO NASA-LIEE 
FORMAT. THIS PROVIDES INPUT TO ‘P‘ OR “P3" 

PROGRAM ACCEPTS MSGS FROM ALL DCP'S AMO 5CPEEHS OUT 
THOSE UHICH ARE NOT NED'S 
OIMEMSION IA(6),lBU2i 


F, 13 ) 

■ 5 . 12 ) 

•5,13) 


CALL APPENDfl5,‘TL0<i*.3,JEft» 5,12> 

IFtlER.NE.DTVPE "OLERMER .. . 

CALL OPEN (S, "USER tSDFM. ICR 1 35 

CALL DFILUCLSaDATSIER) 

CALL CFILU(*LsaDAT*,2,lER) 

CALL 0PeN(e,'LS2DAT*,3.1£R • . 

CALL APPEND (7,"ST0RAQE'’,3,IER> 10O‘ 


IFdER.NE.DTVPE “OPENERR' 

WPITEl-ie.tO) 

URITE( 15,10) 120' 

10 FORPtATi/M 

lCOUNT-0 

ICT*0 • 110 

IYR»6 

' READ BINARY (5)IN0NTH, IDAY,IHR,iniM,ISEC ;GET STARTING TIME 
ISEC-ISEC+-601IMIM 

1 READ BIHARV(5)IA,1S0 ;GET DATA AND SECONDS COUNTER 

ISC- ISCvISEC-ICTtf 3600 R 

IF(I5C.LT.3600)'iO TO 4 

ICT-IC.Tl-1 

ISC*I5C-ICTt3600 

rUDm rUP 4 - 1 

ISr.-I 5 C- 6 LUIpnN 

IFanrtMT*6<5iGO to 

IHP-IHR+1 

iraHF.LT. 24)00 TO 5 

IIiAV-IDAV+l 

tHR «0 

5 CONTINUE 

iF(iAu;.EO.O)oo TO loe 
ICOUNT-ICOUNT+l 


DO ae i-3>6. 
ID»Ittrn,ftND.377IC 
IC»ISHFTriRiI),- 8 » 
IC-ICtrtND,377K 
J-I.ta -1 
IEU)-IC 
j* j+i 
IB<J)-ID 
30 CONTINUE 


ICHK*IA(S >.AMD.20O0«K 

IA<2)-XA{2KAND.7777K jSTfilP EVERTHIHG BUT THE DCP * 

IFaA<2).EO.6063K.OR.lA(2>.EO.65i34K)GO TO 25 
IFaA(2).LT.70O«iaGO TO I 
IF(IA12).E0.7627K)C0 TO 1. 

IFaA(2).Ea.75l4K)G0 TO 1 
IF(IA(2).EQ;7346IOGO TO 1 

8S IER-822^ jASSUnS WO ERROR, OUTPUT A BUNK 


IFciCHK.GT,l)I£R- 17696 *; IF ERR IS FLA(iGED OUTPUT rtN ‘E-^ 
IJRXTEf 6,35 JlYF,inONTH^IDAY,:HR„iniNp ISC, lER,IA(2),vIB( 

UPtTE( 10 , 35 )XVR,r^l 6 NTH,IDf^V,IHP,IMlN,ISC,IER.rA^a),CIB^IC).K 

URlTEaS,35)IVP,inONTH*IDHY,IHR,imN,ISC,IER,IACa),(lB(K),K 

URlTE*'7,3S)IYR,IM0MTH,ID»4Y,lHR,It1ir<,ISC,lER,IA(2), (IBCIO,K- 

CONTINUE 

FORnATii>;>ii, 5 iE,Ai,oi 4 ,eon) 

GO TO I 


CALL CLOSED 5, lEP) 

CALL CL05E*6,LER) 

URITEaO, 150 )I COUNT 
WRITE (15, 120 JlCOUHT 

FORMAT (- TOThL HUMPER OF MESSAGES - M3) 

URlTEa 0 , 110 j 

URITE<15,110j 

FORMAT (/ V/ » 

CALL CHAINS ‘‘P3.3UMEP' 

IFdER^NEanVPE -ODERMEP 

CALL EXIT 

end 
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miGlNAlI T AGE IS 

iOF JOOK QUALEK9 - E-1 5 


Ttt 40 

C TDB -I DEC 75 PROGPr<n TO CONUERT DECINrtL AZIPIUTH AND . 

C ELEI'hTION angles to BCD SUITABLE FOR' INPUT TO 
C ■ INTERFhCE 50 

c TO 1S4S digital couparator 

COKPILER DOUBLE PRECISION 

DIMENSION Kir8),X8<8>,IBC0t«3)»IBCDL(a) XOO 

COMMON IBLLVBCDHIQHt 10 >>BCDL0U(B)«IST0U(5> 

data BcrHiGH/aeo.,ioe..80.,40.,ao.,ie.-,8.,4.*a.«i./ 

•: ECD TABLE 

DAT A BCDLOU*^ «Sf *4^ »3> *1# *08# • 04 # *03# *01^ 

C STOU ANGLES AND FILE ENDER 261 

DATfl ISTOU/'O,O,a0OK,eilK».'-l/ 


CALL 0PENC5,*PTAEM#IER) 
CALL 0PENC6, ’BCDAZEL'^a^IER) 
CALL 0PEN(7, 'HAHGLES^a.IER) 
IFaER.NE.DTVPE *0£MER 

ru«-i 

Dl«l 


R€AD(5.6,END»500)TSIMCE ,READ 1ST REC ON FILE 

READ <5,7)DT,Xa 

F0RMATIF14.9) 

F0RMAT^F5.1/-2ax,F6.a)) 

IFIM»DT 

'JRITE16,11)TS1NCE,DI 

ICT>1 

IFLO«0 

FORMAT' F14.9,F6.1) 

GO TO 109 

IFLG-e 

IFfABStXSa >-Xl( 1 ) ).GT.390>IFLG«1 
IF<IFLG.EQ.i)X8fU-Xa(l)-360. 

■XIIHC* ( Xa< 1 )-Xl ( 1 ) )/DT 

xaiHc-(xa(a>-xi(2))/DT 

DO 105 INC*1,IFIH 
DO 100 1-1,2 lAZ ,EL 
X-XKI; 

IHIOH-0 

ILOU-0 

DO 30 J-1,10 

-n‘lx®t!iGH<J).LE.X>g) TO 20 , CHECK BCD TABLE 

^.0 TO J0 • SMALLER THAN UALUE IN TABLE 

X*:;-ECDHIGHCJ; ;>-OALUE IN TABLE 
II»C+8-(J'lue 


C PRECIOUS LINE IS NAPPING FRON DO LOOP INDEX TO BIT POS 
Oft order on p.9-11 of for 


CONTINUE 


FORT IM 


DO EO J-1,S 

IFtSCDLOUU>.LE.X^GO TO 40 
GO TO 50 


• X-X-BCDLOU(J) 

II-J+6-(J-lilc2 ;BIT MAPPING 
CALL ISET(ILOU,lt). 

CONTINUE 

ISCDHCn-lHlCH 

IBCDL(I)>ILOU 

CONTINUE 

WRITE BINARVC6)IBCDLa),IBCDHCl J,IBCDL<2),lBCDHf2 ' 

ICT-ICT+1 

X1(1)<X1(1HX1INC 

IFCXi(n.GE*e>GO TO 201 

X1(1)-X1(1H360. 

Xl(2)<Xl(2)f'XaiNC 

IF(X1(1)«GT.360. )Xl(l)»XlU)-360. ;PROB HOT NEEDED, B 


CONTINUE 

IF( IFIG.EO. 1 )X3( 1 )-Xaa H360, 

XKi)-xa(i) 

xi(2)-xa(a> 

R£AD(5.110)X2 

F0RNAT(a(lX#F6.a>) 

!F(Xaa).EQ.0,AND.X2(a).E0.9O)G0 TO 120 
GO TO 15 

WRITE BlNARy <6)IST0U 
WRITE<7,700)ICT 
FORMAT ( IX, IG) 

CALI CL0SE(5,IER) 

CALL CL0SE(6,1ER) 

CALL CL0SE(7,IER) 

CALL CHAIN (« LSI. SUMER) 

IF< lER.NE.l )TVPE'INTERP1ERR‘ , lER 

STOP 

END 



APPENDIX F - LITERATURE CITED AND RELATED DOCUMENTS 


1. Escobal, P.R.* "Methods of Orbit Determination"'. 

New York: John Wiley and Sons, 1965 

2. ‘ "How to Use the NOVA Computers". Data General Corporation 

Southboro, Massachusetts 

3. "Real Time Disc Operating System", Revision 3 or higher. 

1 

Data General Corporation, Southboro, Massachusetts 

4. "Fortran IV", Data -General Corporation, Southboro, Massachusetts 


*Ref erred to in text. 
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THE LANDSAT SATELLITE 

and 

FLOOD CONTROL IN NEW ENGLAND 

JUNE 1976 


HISTORY AND BACKGROUND 


Since the Industrial Revolution in the 1800's, the rivers 
of New England have been developed to supply water for 
power and transportation. As new means of transportation 
became more economical both railroad and highway systems 
were built along the banks of the rivers to service the ex- 
panding needs of the industrial, commercial and urban cen- 
ters. Structures, such as buildings, roads, bridges and 
dams have restricted floodways to such an extent that con- 
siderable property and environmental damages have occurred 
during moderate and major floods. Notable floods of Novem- 
ber 1927, March 1936, September 1938 and August 1955 have 
demonstrated the need for flood control to prevent these 
natural catastrophes. 
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At the direction of Congress, the U.S. Army Corps of Engi- 
neers developed a comprehensive plan of protection for each 
river basin after a careful analysis of all water resources. 
Protective works generally consist of a combination of chan- 
nel improvements, dikes and/or floodwalls at major damage 
centers augmented by upstream, flood control reservoirs. 
Many of these reservoirs contain additional storage reserved 
for other uses such as water supply, conservation and rec- 
reation. The Corps has built 35 flood control reservoirs, 
37 local protection projects and 4 hurricane barriers in 
New England at a total investment of over $350 million. 



BALL MOUNTAIN DAM AND RESERVOIR JAMAICA, VERMONT 

To achieve optimum operating benefits from this comprehen- 
sive protection system, the New England Division requires 
hydrologic data such as river, reservoir and tidal levels, 
wind velocity and direction, barometric pressure and pre- 
cipitation. 

In the past this data was collected from field observation 
and relayed via telephone or voice radio. It took several 
hours to compile and assess the data in this manner. With 
the need for timely and reliable information increasing, the 
Corps began development of new methods of data collection. 

n- F - H 


LANDSAT 


Since July 1972, LANDSAT has been relaying river stage, 
precipitation and water quality data from DCP's via the 
Goddard Space Flight Center to the U.S. Army Corps of 
Engineers, New England Division, in near real time. This 

is the first resources satellite designed to obtain data 
from the planet Earth exclusively for planning, design, 
operations and research of land and water resources. 


THE NED GROUND RECEIVE STATION 

Since any operational satellite configuration should include 
ground receiving stations at all major user locales, NED, 
with NASA support has constructed and is now operating an 
inexpensive semiautomatic and easily maintained ground re- 
ceive station as a follow-up to its original study. The 
Division is now able to receive hydrometeorological data 
from data collection platforms in the field directly at its 
headquarters in Waltham, Massachusetts with no time delays. 
The software to drive the antenna system has been developed 
with the intention that the antenna operate in an unattended 
mode automatically over nights and during weekends and hol- 
idays, with a computer controlling all processes. A dia- 
gram of the overall facility is shown. 



SED GROUND RECEIVING STATION DIAGRAM 


In 1970, the Automatic Hydrologic Radio Reporting Network 
was placed in operation. This ground-based radio relay 
system consists of 41 remote reporting stations, and a cen- 
tral control at Division Headquarters in Waltham, Massachu- 
setts. This network, under computer programmed control, 
collects and analyzes, in real time mode, information which 
is essential for flood regulation. The remote reporting 
stations are strategically located in five major river ba- 
sins and at key coastal points, with each contributing to a 
detailed, comprehensive hydrologic picture. 

In June 1972, NASA entered into a contract with the Corps 
for an experiment to study the feasibility of using the 
Earth Resources Technology Satellite (ERTS or LANDSAT) for 
collecting environmental data from Data Collection Plat- 
forms (DCP's) which are installed at 27 locations through- 
out New England. Many are situated at existing U.S. Geo- 
logical Survey gaging stations. 
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DEPARTMENT OF THE ARMY 
NEW ENGLAND DIVISION, CORPS OF ENGINEERS 
WALTHAM, MASS. 


A major objective of the program has been to compare the 
cost, reliability, and operational effectiveness of the 
LANDSAT Data Collection System with the existing NED radio 
network. 

Data collection platforms tested by the Corps have per- 
formed successfully in all seasons including the winter 
months and also during significant flood events, transmit- 
ting near real time operationally useful data for our flood 
fiahtina missions. 


The satellite proved invaluable in April and early May of 
1973 and 1974, monitoring flooding in Maine Rivers. LANDSAT 
\ relayed data from five remote river points in that state to 
aid the New England Division in the coordination of the 
flood emergencies. 


The successful testing of the LANDSAT Data Collection Sys- 
tem at the New England Division should encourage serious 
consideration of the institution of an operational satellite 
data relay system on a Corps -wide basis. System analysis 
is being performed to refine cost data and to articulate the 
data col lection needs of Corps users. 



TRACKING ANTENNA 
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DATA COLLECTION PLATFORM 
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The New England Division is also studying imagery regularly 
collected by LANDSAT to determine the usefulness in plan- 
ning, designing, and managing water resource systems. It is 
well established that such imagery is suited to measuring 
areal extent of ice, snow, and open water, and for estimat- 
ing moisture regimes. Our studies involve computer analysis 
of scenes and will explore indirect methods of calculating 
other hydrologic parameters as well. 



IMAGERY PHOTO TAKEN FROM LANDSAT 


FLOOD CONTROL OPERATIONS 

Data received at the New England Division's Reservoir Con- 
trol Center from either the Automatic Hydrologic Radio Re- 
porting Network or the LANDSAT Data Collection System is 
compiled by computer. This is augmented by information from 
other sources such as the National Weather Service Meteo- 
rologic and River Forecast Offices and the U.S. Geological 
Survey, Experienced engineers and hydrologists at the 
Reservoir Control Center analyze the data for timely opera- 
tion of dams and hurricane barriers, and then issue instruc- 
tions to operating field personnel. 

Flood control reservoirs, local protection projects and 
hurricane barriers built by the Corps in New England have 
been responsible for prevention of about $300 million in 
flood and storm damage. 
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